> 



Structure of the axial-vector meson D s i(2460) and the 

STRONG COUPLING CONSTANT Qd s1 D*K WITH THE LIGHT-CONE 

QCD SUM RULES 

Z. G. Wang0 

Department of Physics, North China Electric Power University, Baoding 071003, 

P. R. China 

iS : 

^ \ Abstract 

In this article, we take the point of view that the charmed axial-vector me- 
. son D s i(2460) is the conventional cs meson and calculate the strong coupling 

[t | \ constant gD sl D*K in the framework of the light-cone QCD sum rules approach. 

The numerical values of strong coupling constants gD sl D*K and gD a0 DK are 
very large, and support the hadronic dressing mechanism. Just like the scalar 
mesons /o(980) and ao(980), the scalar meson D s o(2317) and axial- vector me- 
son L)<ji(2460) may have small cs kernels of the typical cs meson size, the 
t^- ■ strong couplings to the hadronic channels (or the virtual mesons loops) may 

result in smaller masses than the conventional cs mesons in the constituent 
' quark models, and enrich the pure cs states with other components. 

\Q ' 

PACS numbers: 12.38.Lg; 13.25.Jx; 14.40.Cs 
^ \ Key Words: £> sl (2460), light-cone QCD sum rules 

jjn ; 1 Introduction 

The two strange-charmed mesons D s0 (2317) and D sl (2A60) with the spin-parity + 
and 1 + respectively can not be comfortably identified as the quark-antiquark bound 
states in the spectrum of the constituent quark models, they have triggered hot 
debate on their nature, under-structures and whether it is necessary to introduce 
the exotic states [HE]- The masses of the -Dso(23i7) and D s i(2460) are significantly 
lower than the masses of the + and 1 + states respectively from the quark models 
and lattice simulations [3]. The difficulties to identify the D s0 (2317) and D s i(2460) 
states with the conventional cs mesons are rather similar to those appearing in the 
light scalar mesons below lGeV. The light scalar mesons are the subject of an intense 
and continuous controversy in clarifying the hadron spectroscopy [I], the more elusive 
things are the constituent structures of the / (980) and a (980) mesons with almost 
the degenerate masses. The mesons D s0 (2317) and D s i(2460) lie just below the DK 
and D*K threshold respectively, which are analogous to the situation that the scalar 
mesons ao(980) and /o(980) lie just below the KK threshold and couple strongly to 
the nearby channels. The mechanism responsible for the low-mass charmed mesons 
may be the same as the light scalar nonet mesons, the /o(600), /o(980), ao(980) 
and i^o(800) [5J EJ [7J [8] . There have been a lot of explanations for their nature, for 
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example, the conventional cs states [9], [10], two- meson molecular states [TT], four- 
quark states [12], etc. If we take the scalar mesons a (980) and / (980) as four 
quark states with the constituents of scalar diquark-antidiquark sub-structures, the 
masses of the scalar nonet mesons below lGeV can be naturally explained [7J [8]. 

There are other possibilities besides the four-quark state explanations, for ex- 
ample, the scalar mesons a (980), / (980), D s0 (2317) and the axial- vector meson 
-D s i(2460) may have bare P— wave qq and cs kernels with strong coupling to the 
nearby thresholds respectively, the S— wave virtual intermediate hadronic states (or 
the virtual mesons loops) play a crucial role in the composition of those bound states 
(or resonances due to the masses below or above the thresholds). The hadronic dress- 
ing mechanism (or unitarized quark models) takes the point of view that the mesons 
/o(980), a (980), D s0 (2317) and _D s i(2460) have small qq and cs kernels of the typical 
qq and cs mesons size respectively. The strong couplings to the virtual intermediate 
hadronic states (or the virtual mesons loops) may result in smaller masses than the 
conventional scalar qq and cs mesons in the constituent quark models, enrich the 
pure qq and cs states with other components [13], [H] . Those mesons may spend part 
(or most part) of their lifetime as virtual KK, DK and D*K states [SI El EH E]. 
Despite what constituents they may have, we have the fact that they lie just a lit- 
tle below the KK, DK and D*K thresholds respectively, the strong interactions 
with the KK, DK and D*K thresholds will significantly influence their dynamics, 
although the decays D s0 (2317) -> DK and D s i(2460) -> D*K are kinematically 
suppressed. It is interesting to investigate the possibility of the hadronic dressing 
mechanism. 

In our previous work, we take the point of view that the scalar mesons /o(980), 
a (980) and Z} s0 (2317) are the conventional qq and cs state respectively, and calcu- 
late the values of the strong coupling constants gf KK, 9a KK, and Qd s qDK within 
the framework of the light-cone QCD sum rules approach [SI E] • The large values 
of the strong coupling constants support the hadronic dressing mechanism. In this 
article, we take the axial- vector meson D sl (2460) as the conventional cs state, and 
calculate the value of the strong coupling constant Qd s iD*k in the framework of the 
light-cone QCD sum rules approach and study the possibility of the hadronic dress- 
ing mechanism in the axial- vector channel. The light-cone QCD sum rules approach 
carries out the operator product expansion near the light-cone x 2 ~ instead of the 
short distance x ~ while the non-perturbative matrix elements are parameterized 
by the light-cone distribution amplitudes which classified according to their twists 
instead of the vacuum condensates [TSJ HE]. The non-perturbative parameters in 
the light-cone distribution amplitudes are calculated by the conventional QCD sum 
rules and the values are universal [TT] . 

The article is arranged as: in Section 2, we derive the strong coupling constant 
9d s1 d*k within the framework of the light-cone QCD sum rules approach; in Section 
3, the numerical result and discussion; and in Section 4, conclusion. 
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2 Strong coupling constant gr> sX D*K with light-cone 
QCD sum rules 

In the following, we write down the definition for the strong coupling constant 
9d s1 d*k , 

(D sl (p + q)\D*(p)K(q)) = -ig D . lD .KV*a<? = -tM Dsl g DslD * KV * a e a , (1) 

where the e a and r] a are the polarization vectors of the mesons D* and D sl (2460) 
respectively. The mass of the _D s i(2460) Md s1 can serve as an energy scale, we 
factorize the Md s1 from the Qd s1 d*k- We study the strong coupling constant gD al D*K 
with the two-point correlation function IV (p, q), 

lWp,g) = ijd*xe-^{0\T{jV(0)jA + ( x )}\K(p)), (2) 
Jj( x ) = n(x)7 M c(x), (3) 
Jfa) = s(x)7 M 7 5 c(x) , (4) 

where the vector current JY(x) and the axial- vector current J^(x) interpolate the 
vector meson D* and the axial- vector meson D s i(2460) respectively, the external K 
state has four momentum p^ with p 2 = m 2 K . The correlation function IV(j>, q) can 
be decomposed as 

IV (p, q) = iRg^v + ni(p M g„ + p^) H (5) 

due to the Lorentz invariance. 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach [17] . we can insert a complete series of intermediate states with 
the same quantum numbers as the current operators JY(x) and Jf}(x) into the 
correlation function IV (p, q) to obtain the hadronic representation. After isolating 
the ground state contributions from the pole terms of the mesons D sl (2460) and 
D* , we get the following result, 

(Ol-ff(O) \D*(q + p))(D*\DsiK}(D sl (q)\J^(m 
[Ml,-{q + pf][Ml ai - q 2 ] 

i9D. 1 D-KfD-fD*M D .M Dml ^ (6) 



[Mh*-{q+v) 2 \ [M 2 Dsl -q 2 ] 
where the following definitions have been used, 

{0\J%(PW) = fD*M D *6^ 

(0|J»|Z? al > = f D3l M DslV ,, 
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here the fo* and fu al are the weak decay constants of the D* and D sl (2460) respec- 
tively. The vector current JY(x) and axial- vector current Jw~(x) have non- vanishing 
couplings to the scalar meson D and pseudoscalar meson D s , respectively, 

(0|jJ(0)|A)(<z)> = footo, 

(0\J$(0)\D s (q)) = ifafy, 

where the fn and fo s are the weak decay constants. The 111 with the tensor 
structure p^q u + p u q^ receives contribution from the mesons D and D s besides the 
D* and Z) s i(2460), we choose the tensor structure g^ u for analysis to avoid possible 
contaminations from the scalar and pseudoscalar mesons. In Eq.(6), we have not 
shown the contributions from the high resonances and continuum states explicitly 
as they are suppressed due to the double Borel transformation. 

In the following, we briefly outline the operator product expansion for the correla- 
tion function n^ w (p, q) in perturbative QCD theory. The calculations are performed 
at the large space-like momentum regions (q + p) 2 <C and C 0, which corre- 
spond to the small light-cone distance x 2 ~ required by the validity of the operator 
product expansion approach. We write down the propagator of a massive quark in 
the external gluon field in the Fock-Schwinger gauge firstly [18J, 

(o|r{ ft (m)^(x 2 )}|o> = • / 'J^ e -*c»™> 



(27T) 
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! " " l kj ~ I dvg s G^{vx x + (1 - v)x 2 ) 



o 

1 M + m 1 



2 (A; 2 — m 2 ) 2 ^ k 2 — m 2 

where the is the gluonic field strength, and the g s denotes the strong coupling 
constant. Substituting the above c quark propagator and the corresponding K 
meson light-cone distribution amplitudes into the correlation function n^p, q) in 
Eq.(2) and completing the integrals over the variables x and k, finally we obtain the 
analytical result, which is given explicitly in the appendix. 

In calculation, the two-particle and three-particle K meson light-cone distri- 
bution amplitudes have been used [15j [161 EH EH 120], the explicitly expressions 
are given in the appendix. The parameters in the light-cone distribution ampli- 
tudes are scale dependent and can be estimated with the QCD sum rules approach 
[T51 [TBI [TBI [JjJJ ED]- In this article, the energy scale \x is chosen to be [i — lGeV. 

Now we perform the double Borel transformation with respect to the variables 
Q\ = —q 2 and Q\ = — (p + q) 2 for the correlation function IT in Eq.(6), and obtain 
the analytical expression of the invariant function in the hadronic representation, 



-v{xx - X 2 )^u \ , (7) 



here we have not shown the contributions from the high resonances and continuum 
states explicitly for simplicity. In order to match the duality regions below the 
thresholds s and s' for the interpolating currents JY(x) and J^{x) respectively, we 
can express the correlation function II at the level of quark-gluon degrees of freedom 
into the following form, 

n = / dsds'- r — , p<yS \ s 2 * , (9) 

J [s-(q + p) 2 ][s'-q 2 ] KJ 

then perform the double Borel transformation with respect to the variables Q\ and 
Q\ directly However, the analytical expression of the spectral density p(s, s') is 
hard to obtain, we have to resort to some approximations. As the contributions 
from the higher twist terms are suppressed by more powers of m 2_^ +up ^2 , the net 
contributions of the three-particle (quark-antiquark-gluon) twist-3 and twist-4 terms 
are of minor importance, about 20%, the continuum subtractions will not affect the 
results remarkably. The dominating contribution comes from the two-particle twist- 
3 term involving the <p p (u). We preform the same trick as Refs. [T8| |2T] and expand 
the amplitude <p p (u) in terms of polynomials of 1 — u, 



N N 



y-v , / -s — III 

k=0 k=0 



^(tt) = E 6 *( 1 - u )* = E 6 *(73^) ' ( 10 ) 



then introduce the variable s' and the spectral density is obtained. In the decay 
B — > XcoK, the factorizable contribution is zero and the non-factorizable contribu- 
tions from the soft hadronic matrix elements are too small to accommodate the ex- 
perimental data [22], the contributions of the three-particle (quark-antiquark-gluon) 
distribution amplitudes of the mesons are always of minor importance comparing 
with the two-particle (quark-antiquark) distribution amplitudes in the light-cone 
QCD sum rules. In our previous work, we study the four form-factors fi(Q 2 ), 
fiiQ 2 )-, 9i{Q 2 ) and giiQ 2 ) of the £ — > n in the framework of the light-cone QCD 
sum rules approach up to twist-6 three-quark light-cone distribution amplitudes and 
obtain satisfactory results [23J. In the light-cone QCD sum rules, we can neglect the 
contributions from the valence gluons and make relatively rough estimations. 

After straightforward calculations, we obtain the final expression of the double 
Borel transformed correlation function U(M 2 ,M 2 ) at the level of quark-gluon de- 
grees of freedom. The masses of the charmed mesons are Md s1 = 2A6GeV and 
Mb* = 2.01GeV, M — ~ 0.45, there exists an overlapping working window for 

the two Borel parameters M\ and Mf , it's convenient to take the value M\ = Mf . 
We introduce the threshold parameter sq and make the simple replacement, 



m ?+ u 0( 1 — "0) m |f m e+"o( 1 ~ u o) ro jc s l) 

IT 1 — > P W 2 — ft _ M' 2 



to subtract the contributions from the high resonances and continuum states 
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finally we obtain the sum rule for the strong coupling constant c/d s \D* k i 



9d s1 d*k 



where 



f D *f Dsl M D *M Dsl 
fx 



cxp 



M 2 Dsl Ml, 



+ 



Mf M 2 



exp 



BB 

Up 



m c m 2 K M 2 . . m 2 K (M 2 + m 2 ) d .. , M 4 d 



m u + m s 
BB 



cxp 



+m 2 K 



M 2 



da* 



PUO 

f K m 2 c m 2 K / dtB(t) 
Jo 



2 du 



d 



l—Cis 



uo-a a 
Uq 



da„ 



+f K m 2 K M z — / da a 
du Jo 



(uof K m 2 K § + fsKm c (j) 3K )(l -a s - a g , a s , a 6 

a g 

1 ~ a - (A\\-V\\)(l-a a -a g ,a a ,ag) 
aa n 



d 



du J 

l—UQ 



UO—Cts 

1— a s 



2a n 



Jk^k^J- j — / da s / da g A\\ (1 — a s — a g , a 8 , a 



a s + a„ — no 



no— "s 



+fKlTT>K / / / da + 



L«9 





2m 2 
d M 2 

d 



UQ—CLg 

2 



da n 



l—UQ J UQ—OLg 



9 J 

1-Q„ 



a 



da= / da 



M 2 



a; 



1 — a — a g , a, a g ) 



-fKm 4 K u 



l—UQ 



du y 

<3>(1 — a — a g , a, a 9 ) 
a„ 



da„ 



UQ-Otg 



da s / da + 



da„ 



l-u Ju —a g 



l-a a 



da Q / da 



/j<-m K / da 9 / dp 

Jl-U 

Am 2 (1 - m ) 2 



1-/3 



da 



$(1 - a - /?, a,/?) 



l-it 



of 



2m 2 c 



M 2 a\ 



(A,, +A ± )(l-a-/3,a,/3) 



d 



diio A 



-«o Jo 



+/xm^- — / da 5 / dp da$(l — a — /?, a, /3) 



1-/3 



u (l - «o) 



Of 



'111 



SB 

M 2 



m 2 + w (l — Wo)^ 
M 2 



2 

K ! 



M 2 + M| ' 
M 2 M 2 



M 2 + M| ' (12) 

here we write down only the analytical result without the technical details B The 
term proportional to the M 4 ^0#-(uo) in Eq.(ll) depends heavily on the asymmetry 



In this footnote, we present some technical details necessary in performing the Borel trans- 
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coefficient ai(/i) of the twist-2 light-cone distribution amplitude <Pk{u) in the limit 
u = |, if we take the value ai(/x) = 0.06 ± 0.03 [19J, [20], no stable sum rules can 
be obtained, the value of the gD sl D*K changes significantly with the variation of 
the Borel parameter M 2 . In this article, we take the assumption that the u and s 
quarks have symmetric momentum distributions and neglect the coefficient ai(/i). 
The existence of such a term is not a bad thing. The -D s i(2460) lies below the D*K 
threshold, it is impossible to measure the strong coupling constant Qd s1 d*k directly, 
the corresponding beauty doublet (0 + , 1 + ) B s0 and B s i may lie above the BK and 
B*K thresholds respectively. Once the experimental data of the B s0 , B s i and the 
related strong coupling constants are available, we can compare the values of the 
9b s1 b*k from the QCD sum rules with the ones from the experiment, and verify 
whether or not the coefficient ai(/i) can be safely neglected. That will put severe 
constraint on the value of the ai(/i). With the simple replacement 



m c - 


-> m b , 


M Dsl - 


"> M Bsl 


Mb. - 


-> Mb* 


fo sl - 




Id* ~ 


- /a- 



in Eq.(ll), we can obtain the QCD sum rule for the strong coupling constant Qb s1 b*k- 



3 Numerical result and discussion 

The parameters are taken as m s = (140 ± 10)MeV, m q = (5.6 ± 1.6)MeV, m c = 
(1.25±0.10)GeV, A 3 = 1.6 ±0.4, f 3K = (0.45 ±0.15) x lO~ 2 GeV 2 , u 3 = -1.2 ±0.7, 
774 = 0.6 ± 0.2, lu 4 = 0.2 ± 0.1, a 2 = 0.25 ± 0.15 [151 IIS H3 US 120], !k = 0.160Ge1/, 



formation which are not familiar to the novices, 



dv / da, 



o 

l 



da s f{v,a s ,a g )—e-xv 



m 2 c + u(l — u)m 



K 



M 2 



S(u - Uo)|«=a a +(l— v)a g 



du I dv I da 



da s f{v, a s ,a g )5 [u - a s - (1 - v)a g 



— exp 

du 



m 2 . + u(l — u)m 



K 



M 2 



5{u - u Q ) 



du / da & 



f(v,a s ,a g ) d 
da„ — — exp 



du 



m 2 + u(l — u)m 2 K 



M 2 



du exp 



m 2 + u(l — v)m 



K 



M 2 



S(u - wo) j- / da s 



da, 



S(u - u ) 
f(v,a s ,a g ) 



exp 



m 2 + ito (1 — uo)m 



M 2 



d 

du j 



da. 



da, 



f(v,a s ,a g ) 
a a 



'UQ—a s ^9 

where the f(v,a s ,a g ) stands for the three-particle light-cone distribution amplitudes. 
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m K = 49SMeV, M Dal (2460) = 2A6GeV, M D * = 2MGeV, f D * = (0.24 ± 0.02)GeV 
HB], and f Dal = (0.225 ± 0.020)GeV [TO]. The duality threshold s in Eq.(ll) 
is taken as Sq = (6.8 — 7.2)GeV 2 to avoid possible contaminations from the high 
resonances and continuum states, it is reasonable for the narrow D s i(2460) = 
(2.6 — 2.7)GeV > Md s1 , furthermore, in this region, the numerical result is not 
sensitive to the threshold parameter s - The Borel parameters are chosen as Mf = 
Mf = (6 — 14)GeV 2 and M 2 = (3 — 7)GeV 2 , in those regions, the value of the strong 
coupling constant Qd s iD*k is rather stable from the sum rule in Eq.(ll) with the 
simple subtraction, which is shown in Figs. (1-2). 

The uncertainties of the four parameters m u , m c , A3 and U3 can only result in 
small uncertainties for the numerical values. The main uncertainties come from 
the seven parameters f 3 x, m s, Id*, Id s1 i ^2, Vi an d ^4, the variations of those 
parameters (except for the 02) can lead to large changes for the numerical values, 
about (5 — 10)%, which are shown in the Fig.l. The uncertainties of the three 
hadronic parameters fox, /d*, /d s i and the two light-cone distribution parameters 
774, can be pined down with the improved QCD sum rules or more experimental 
data, however, it is a difficult work. 

Taking into account all the uncertainties from the thirteen parameters m s , m q , 
m c A3, fen, <^3, T]i, tt2, /d* , /d s1 , so and M 2 , finally we obtain the numerical 
result of the strong coupling constant, 

9d s1 d*k = (10.5±3.5)Gey, (14) 
9d. iD .k = 4.3 ±1.4, (15) 

the uncertainty is large, about 30%. The large values of the strong coupling constants 
9d s0 dk (QDsoDK = (9.6 ± 2A)GeV [B]) and Qd s \D*k obviously support the hadronic 
dressing mechanism 0, the scalar meson D s q(2?>17) and axial-vector meson D s i(2460) 
(just like the scalar mesons /o(980) and a (980), see Ref. [5]) can be taken as hav- 
ing small scalar and axial-vector cs kernels of typical meson size with large virtual 
S'-wave DK and D*K cloud respectively. In Ref. [24] . the authors analyze the unita- 
rized two-meson scattering amplitudes from the heavy-light chiral Lagrangian, and 

3 Here we will take a short discussion about the hadronic dressing mechanism |13[ I14j. one can 
consult the original literatures for the details. In the conventional constituent quark models, the 
mesons are taken as quark-antiquark bound states. The spectrum can be obtained by solving 
the corresponding Schrodinger's or Dirac's equations with the phenomenological potential which 
trying to incorporate the observed properties of the strong interactions, such as the asymptotic 
freedom and confinement. The solutions can be referred as confinement bound states or bare 
quark-antiquark states (or kernels). If we switch on the hadronic interactions between the confine- 
ment bound states and the free ordinary two-meson states, the situation becomes more complex. 
With the increasing hadronic coupling constants, the contributions from the hadronic loops of the 
intermediate mesons become larger and the bare quark-antiquark states can be distorted greatly. 
There may be double poles or several poles in the scattering amplitudes with the same quantum 
number as the bare quark-antiquark kernels; some ones stem from the bare quark-antiquark kernels 
while the others originate from the continuum states. The strong coupling may enrich the bare 
quark-antiquark states with other components, for example, the virtual mesons pairs, and spend 
part (or most part) of their lifetime as virtual mesons pairs. 
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f =0.45*1 0" 2 GeV 2 ; 
f 3K =0.60*10" 2 GeV 2 ; 
f 3K =0.30*10" 2 GeV 2 . 
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Figure 1: The gD sl D*K with the parameters f 3 x, m s> Jd*> Jd s1 , <^4 and 774 respec- 
tively. 
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Figure 2: The Qd s \D*k with the Borel parameter M 2 . 



observe that the scalar meson D s o(2317) and axial-vector meson Z) s i(2460) appear 
as the bound state poles with the strong coupling constants Qd s oDK = 10.203GeF 
and Qd s \D*k = 10.762Gey. Our numerical results gD s0 DK = (9.6 ± 2A)GeV and 
9d s1 d*k = {!0.5±3.5)GeV are certainly reasonable and can make robust predictions. 
However, we take the point of view that the meson D s o(2317) (D s i(2460)) be bound 
state in the sense that it appears below the DK (D*K) threshold, its constituents 
may be the bare cs state, the virtual DK (D*K) pair and their mixing, rather than 
the DK (D*K) bound state. In Ref . }25] . the authors take the point of view that 
the .D s o(2317) is the scalar cs meson and calculate the mass Md s0 with the QCD 
sum rules approach by taking into account the contribution of the DK continuum, 
the effects of the DK continuum can pull the mass down remarkably, and the value 
of the Mo s0 is in good agreement with experimental data. Our numerical values 
of the strong coupling constants Qd s $dk and Qd s1 d*k are approximately equal, the 
spin symmetry of the heavy quarks works rather well, the contribution of the D*K 
continuum may pull the mass Md s1 down remarkably. One can analyze the value 
of the Md s1 in the framework of the QCD sum rules with the axial-vector current 
Jn{ x ) by including the contribution of the D*K continuum. 
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4 Conclusion 



In this article, we take the point of view that the charmed mesons D s0 (2317) and 
.D s i(2460) are the conventional cs mesons and calculate the strong coupling constant 
9d s1 d*k within the framework of the light-cone QCD sum rules approach. The nu- 
merical values of the strong coupling constants Qd s \D*k and Qd s qDK are compatible 
with the existing estimations, the large values support the hadronic dressing mecha- 
nism. The uncertainty of the value of the Qd s1 d*k is large, about 30%, it comes from 
the uncertainties of the thirteen parameters m s , m q , m c , A3, fax, ^3, fyt, ^4, 0-2, I'd*, 
Id s1 , So and M 2 , while the main uncertainty comes from the seven parameters f^K, 
m s, Id*, /d s i, a 2, and UJ4, the variations of those parameters (except for the a 2 ) 
can lead to large changes for the numerical values, about (5 — 10)%, refining those 
parameters is of great importance, improved QCD sum rules and more experimental 
data may pin down the uncertainties. 

Just like the scalar mesons /o(980) and ao(980), the scalar meson D s o(2317) and 
the axial- vector meson D sl (2460) may have small cs kernels of typical cs meson size. 
The strong couplings to virtual intermediate hadronic states (or the virtual mesons 
loops) can result in smaller masses than the conventional + and 1 + mesons in the 
constituent quark models, enrich the pure cs states with other components. The 
-D s o(2317) and D s i(2460) may spend part (or most part) of their lifetimes as virtual 
DK and D*K states, respectively. 
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Appendix 

The analytical expression of the II at the level of the quark-gluon degrees of freedom, 



m u + m s Jq AA Jq Jq AA 2 

A(u)m 2 K d 



+ lfJ\ u { M u)^o g AA + 4 du 
+m 2 K f dv f da„ / da s 



1 ml 
AA + AA 2 



'0 Jo JO 



~^j^2 \u=a s +{l-v)a g 

— - dv da g I da s [(1 - 2v)A\\ - V\\] (1 - a s - a g , a s , a g ) 

A Jo Jo Jo 



d 1 

\u=a s + (l— v)a g 

■1 r l rl-a 



du AA 



+fK m K dv da g / da s / da&(l 
Jo Jo Jo Jo 

J 4 Am 2 c d 1 1 
\ AA 2 ~ AA* +U d^AA 2 ~\ „ , 

V j u=a a + (l— v)a 



a 

'fKm-K / dvv / da g / df3 I da&(l — a — (3, a, (3) 
Jo Jo Jo Jo 



J 4 Am 2 c d 1 1 
{AA 2 ~ AA* + U \Tu~AA 2 ) u=1 _ 



+%fKrn 2 c m 4 K I dvv / da g / da s / da{A\\ + A±)(l — a — a g , a, a g ) 
Jo Jo Jo Jo 



1 

u=a s +(l-v)a g 



AA 3 



a — (3, a, (3) 



AA 3 
where 



r 1 r 1 r a g r^-P 

SfxTnlm^ / dvv 2 / da g / d(3 I da(A\\ + A±)(l 
Jo Jo Jo Jo 

u=l—va„ 1 (16) 



AA = ml — (q + up) 2 , 
$ = A\\ + A ± -V\\-V ± . (17) 
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The light-cone distribution amplitudes of the K meson are defined by 



\u(0)w 5 s(x)\K(p)) = if K Pn I due~ tup ' x <{ <f> K (u) + —^-A(u) 



o 



2 2 

m K x 



+fKm 2 K7 p^- I due-^ x B{u) 
2p ■ x Jo 



(0\u(0)ilBs(x)\K(p)) = fK ™ K [ due-**-*<p p (u) , 

m s + m u Jo 

(0\u(0)a^'j 5 s(x)\K(p)) = i(p M a;„ - p v x^) f R ™ K I due~ mp ' x (f a (u) , 

o{m s + m u ) Jo 

(0\u(0)a a ^ 5 g s G flu (vx)s(x)\K(p)} = f 3K {{p^p a g^ - p v p a 9^p) ~ (P^PpgL 

-PuPf>g}a)} J Va^K^e-^^K 
(0\u(0)>y fl >y 5 g s G a p(vx)s(x)\K(p)} = p PaXp - V<iXa f K m 2 K 

p ' X 

+ fKm 2 K (p p g a ^ - p a gpn) 
J Da i 4(« i )e" ip " (as+TO9) , 

(0\u{0)~f^g s G al3 {vx)s{x)\K(p)) = p PaXl3 - VpXa f K m 2 K 

J Va i V ll (a i )e- ip < as+va ^ 
+fKm 2 K (pf 3 g a ^ - p a gp,j) 

J Da^ifa.je-^^^' , (18) 

where the operator G a p is the dual of the G a/ 3, G a/ 3 = \£ a pnuG^ v and T>a.i is defined 
as Va-i = da.idot2dotz5{l — «i — a 2 — 0:3). The light-cone distribution amplitudes are 
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parameterized 

<j) K {u,n) = 
<P P (u,n) = 



as 



(p a (u,[l) 
<l>3K(<Xi, fJ>) 

V\\(ai,n) 

A\\(Oi,fi) 

V ± (aii,n) 

A ± (ai,fi) 
A(u,n) 



9k(u,h) 
B(u,n) 



6u{l - u) |l + a 1 c}(2u - 1) + a 2 c\{2u - 1) + a A Cl(2u - 1)} , 

i + jso^-^jcf^-i) 



6u(l-u) < 1 + 



1 7 2 3 2 ' 

1 



C?(2u-1) 



360a u a s a fl jl + A 3 (a u - a s ) + cj 3 -(7« 9 - 3) 

120a u a s a 9 (w 00 + Vi (3ai g - 1)) , 

120a u a s a g a w (a s - a u ) , 

-30a 2 {h 00 (l - a g ) + h 01 [a g (l - a g ) - Qa u a s 



+h 



10 



20 

a 2 + 2O773 + —774 
y 



+ 
+ 



(a« - 


Ota) 




— u) < 


f 16 

1 15 


24 
+ 35 


1 

-I5 + 


1 


7 


16 " 


"27 f ' 
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4 



10 
27 



C*(2u-1) 



210 



a 2 



135 



%^3 



5 



a 2 + 21774(^4 



C 4 3 (2«-l)| + 

{2m 3 (10 - 15m + 6m 2 ) logu + 2w 3 (10 - 15u + 6u 2 ) logw 
+tm(2 + 13mm)} , 

l + 52 C|(2M-l)+ 5 r 4 C|(2M-l), 
9k(u,h) - <M«,«) , 



(19) 
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where 



7 Vi 

n oo — v oo — — y ) 

21 9 
aio = y^W4 - — a 2 , 

21 

^10 = -^V^4 , 
o 

fc 7 3 

«oi = - — a 2 , 

7 3 

18 20 
#2 = 1 + y a2 + 60% + y 7/4 , 

9 

^4 = -^g«2 - 6773^3 , (20) 

- - - f + 2 

here C 2 2 , C 4 2 and C 2 2 are Gegenbauer polynomials, 773 = ^ m ^ 2 ms and p 2 = -^f- 

nsi hhi hhi eq]. 
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